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Stefano Varricchio suddenly passed away on August 20th 2008, shortly after 
this paper has been completed. We will remember Stefano as the best friend 
of us and as an outstanding researcher. Working with Stefano was always an 
enthusiastic experience both for his beautiful and original ideas and for his 
scientific rigueur. 

Abstract 

Let X be a subset of N l or Z*. We can associate with X a function 
Qx ■ N — > N which returns, for every (ni, . . . ,rit) £ N , the number 
Qx (ni , . . . , Tit) of all vectors x € X such that, for every i = 1, . . . , t, \xi\ < 
rii. This function is called the growth function of X. The main result of 
this paper is that the growth function of a semi-linear set of N 4 or 1} is 
a box spline. By using this result and some theorems on semi-linear sets, 
we give a new proof of combinatorial flavour of a well-known theorem by 
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Dahmen and Micchelli on the counting function of a system of Diophantine 
linear equations. 

1 Introduction 

We study some counting problems on semi-linear sets of the free commutative 
monoid N 4 and of the free commutative group Z* . The notion of semi- linear set 
has been widely investigated in the past because it plays an important role in 
the study of several problems of Mathematics and Computer Science. In this 
context, it is worth to mention a well-known theorem by S. Ginsburg and E. 
H. Spanier [10] (see also [9] [13]) that establishes a strong connection between 
semi-linear sets, rational sets on N and on Z and Presburger definable sets. In 
[S] , Eilenberg and Schutzenberger extended the connection between semi-linear 
sets and rational sets to every finitely generated commutative monoid. 

It is now convenient to introduce two concepts that are central in our work: 
those of quasi-polynomial and box-spline. A quasi-polynomial is a map F : 
N* — > N, defined by a finite family of polynomials in t variables x%, . . . ,Xf, with 
rational coefficients: 

{P(d u d 2 ,-,d t ) I di,...,dt €N, 0<di <d}, 

such that every polynomial of the family is indexed by a vector (d\, d,2, ■ • ■ , dt) 
whose components are remainders of a fixed positive integer d. Then, for every 
(m, . . . ,n t ) <E N*, the value of F computed at (m, . . . , n t ) is given by: 

F(m, ...,n t )= P( dl ,d 2 ,---,d t )(ni, n t ), 

where, for every i = 1, . . . , t, dj is the remainder of the division of Uj by d. 

In this paper, we shall use a particular kind of box-spline specified as a 
map F : N* — > N such that there exist a partition of N* into a finite num- 
ber of polyhedral conic regions R\, . . . , R s , determined by hyperplanes, through 
the origin, with rational equations, and a finite number of quasi-polynomials 
pi, . . .p s , where, for any (m, . . . , rit) one has: 

F(m, ...,n t )= Pj(ni, . . . ,n t ), 
where j is such that (rii, . . . , n t ) e Rj. 

In the sequel, we shall use the notion of box-spline in the sense specified above. 
One interesting application of these notions was first given by E. T. Bell in [TJ, 
by proving that the counting map of a diophantine linear equation is a quasi- 
polynomial. Later, this result was extended in a paper by Dahmen and Micchelli 
[5] , where it has been proved that the counting function of a Diophantine system 
of linear equations can be described by a set of quasi-polynomials, under suitable 
conditions on the matrix of the system. Recently this result has been object of 
further investigations in [6l [TJ [15] where important theorems on the algebraic 
and combinatorial structure of partition functions have been obtained. In formal 
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language theory, the notion of quasi-polynomial has been used to describe the 
counting and the growth function of regular languages [14] and of bounded 
context-free languages [31 H]. I n particular, in 0], a result perhaps of some 
interest is that the Parikh counting function of a bounded context-free language 
is a box-spline. Moreover, starting from the grammar that generates such a 
language, one can effectively construct a box-spline that describes its Parikh 
counting function. 

In this paper, we introduce and study the notion of growth function of semi- 
linear sets both over N as well as over Z. More precisely, we define the growth 
function of a subset X of N 4 or Z* , t > as the function fx ■ N* — ► N which 
associates with non negative integers n\, . . . ,nt, the number fx(ni, . . . ,n t ) of 
all the elements (rti, . . . , n t ) £ X such that: 

\xi\ < ni 

< 

\x t \ < n t . 

The function fx seems to be a natural way to count the number of points of 
X. Indeed, in case X lies in N* (resp. Z'), it counts the number of points of X 
lying in larger and larger hyper-parallelepipeds, starting from the origin of N* 
(resp. centered at the origin of Z*). In this paper, by using the combinatorial 
techniques developped in [4] , we prove that the growth function of every semi- 
linear set of N* or Z', t > 0, is always a box spline. Furthermore, this box-spline 
can be effectively computed. Moreover, by using such techniques and some deep 
results from the theory of semi-linear sets, we give a new proof of the theorem 
of Dahmcn and Micchelli mentioned above. 

The paper is organized as follows. In Section 2, after introducing some basic 
notions and results, we recall the definition of growth function of a semi-linear 
set. In Section 3, we give a combinatorial proof of the fact that the counting 
function of a system of diophantinc linear equations with coefficients in N is a 
box spline. This will be a crucial tool to prove the results of the subsequent 
sections. We remark that the proof of the result mentioned above has been 
already published in [3] and we reproduce it here for the sake of completeness. 
In Section 4, we study the growth functions of semi-linear sets and, as the main 
result of this section, we prove that such functions are box-spline. Finally, in 
Section 5, we give a proof of combinatorial flavour of the quoted theorem by 
Dahmen and Micchelli. 

2 Preliminaries and basic definitions 

The aim of this section is to recall some results about semi-linear sets of the 
free commutative monoid and the free commutative group. For this purpose, we 
follow [2]. The free abelian monoid and the free abelian group on k generators 
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are respectively identified with N fe and Z fe with the usual additive structure. 
The operation of addition is extended from elements to subsets: if X, Y C N* 1 
(resp. I,rc 1 k ), X + Y CN k (resp. X + Y C Z fc ) is the set of all sum a: + y, 
where x £ X,y £ Y . It might be convenient to consider the elements of N k and 
Z fc as vectors of the Q fc -vector space Q fc . Given v in N fc or in Z fc , the expression 
Nv stands for the subset of all elements nv, where neN. This expression can 
be extended to Zv, whenever v is in Z k . Let B = . . . , b n } be a finite subset 
of Z . Then we denote by B® the submonoid of N k generated by B, that is 

B® = b® -\ lb® = {mjiiH h m„6„ | m, € N}. 

In the sequel, the symbol K stands for N when it concerns the free abelian 
monoid N fc and for Z when it concerns the free abelian group Z fe . The following 
definitions are useful. 

Definition 1 Let X be a subset of Z fc (resp. N k ). Then 

1. X is K- linear if it is of the form 

n 

a + ^Kbi, a, 6 i eZ fe , (resp.N k ) 1 i = l,...,n; 

i=i 

2. X is K-simple if the vectors bi are linearly independent in Q k , 

3. X is K-semi-linear if X is a finite union ofWL-linear sets; 

4. X is semi-simple if X is a finite disjoint union of K-simple sets. 

Remark 1 In the definition of simple set, the vector a and those of the set 
{&!, . . . , b n } shall be called a representation of X. 

There exists a classical and important connection between the concept of semi- 
linear set and the Presburger arithmetic. Denote by Z = (Z; =; <; +; 0; 1) and 
by M = (N; =; +; 0; 1) respectively the standard and the positive Presburger 
arithmetic. Given a subset X of N fe (resp. Z fc ), we say that X is first-order 
definable in M (resp. Z), or a Presburger set of N k (resp. Z fc ), if 

X = {(xi, ■ ■ ■ ,Xk) I P(xi, . . . ,Xk) is true}, 

where P is a Presburger formula (with at most k free variables) over N (resp. 
Z). Ginsburg ans Spanier proved in [10] the following result for N fe that holds 
for Z fc also. 

Theorem 1 (Ginsburg and Spanier, 1966) Given a subset X o/N fe (resp. 1r), 
the following assertions are equivalent: 

1. X is first- order definable in AT (resp. Z); 

2. X is N-semilinear in N k (resp. Z k ); 
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3. X is N-semisimple in N k (resp. 1 k ). 

Remark 2 A Z-semilinear set of Z fc is always N-semilinear in Z fe . 

Remark 3 Theorem [T] is effective. Indeed, one can effectively represent a N- 
semi-linear set X as a semi-simple set. More precisely, one can effectively con- 
struct a finite family {Vi} of finite sets of vectors such that the vectors in Vi 
form a representation of a simple set X{ and X is the disjoint union of the sets 
Xi. 

Remark 4 Given a monoid M, a subset of M is rational if it is obtained from 
finite subsets of M by applying finitely many times the rational operations, that 
is, the set union, the product, and the Kleene closure operator. Obviously, a 
semi-linear set of N fc or 1 k is rational but one can prove that the opposite is 
true (see [13]). Therefore, Conditions 2 and 3 of Theorem [1] and the property of 
rationality are equivalent in N fc or X k . We recall that the previous equivalence 
has been proven in the larger context of finitely generated commutative monoids 
by Eilenberg and Schutzenberger [5]. 

The following definition is important. 

Definition 2 Let X be a subset of "P. We associate with X a function 

Qx ■ N* — ► N 

defined as: for every (m, ...,nt) S N*, Qx(ni, ■■■,nt) is the number of elements 
(xi, Xt) £ X such that: 

\x\\ < n\ 

i ■ ■ (1) 

\x t \ < n t . 

The function Qx is called the growth function of X . 

3 Counting the non-negative solutions of sys- 
tems of Diophantine linear equations with pos- 
itive coefficients. 

As we have said in the introduction of this paper, our main goal is to give an 
exact description of growth functions of semi-linear sets of Z fe . For this purpose, 
now we develop a combinatorial tool to deal with such functions. This tool con- 
cerns systems of Diophantine linear equations. We start by considering systems 
of Diophantine linear equations with positive coefficients. Given a system of this 
kind, we can define a function that returns, for any vector of non homogenuous 
terms, the number of non-negative solutions of the system. We will prove that 
this function is a box-spline. 
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Lemma 1 Let q(x\, . . . ,Xt,x) be a polynomial in t + 1 variables with rational 
coefficients and let 

F : N* x {{-1}UN} — > Q 



be the map defined as: 

Ea=o,..,j) li x U ■■■,x u \) x > 0, 
x=-l. 



F(x-l, ...,x t ,x) 



There exists a polynomial p(x±, . . . , Xt, x) in t + 1 variables with rational coeffi- 
cients such that, for every (x\, . . . , x t , x) G N* x {{—1} U N}, one /ias: 

F{x u ...,Xt,x) =p(xi,...,x t ,x). 

Proof. Write q(x±, . . . , xt, x) as: 

a a + aix + ■ ■ ■ + a n x n , (2) 

where, for every i — 0, . . . , n, at is a suitable polynomial in the variables 
Xi, . . . , Xt with rational coefficients. By Eq. I[2]). if a; > 0, for every Xi, . . . , x% € 
N, one has: 

F(x!,...,Xt,x) = ^ gr(xi,...,a; t ,A) = a J ' X! A ^ ' ( 3 ) 

A=0,...,x ,j'=0,...,n V A=0,...,x / 

On the other hand, by using a standard argument (cf 4J, Appendix), one can 
prove that, for any j £ N, there exists a polynomial Pj{x) with rational coeffi- 
cients in one variable x such that: 

3.1) for any igN, pj(x) = Y,\=o,...,x A ^ 

3.2) Pj -(-l) = 0. 

For any j = 0, . . . , n, let pj be the polynomial defined above and let p — 
p{x\, . . . , Xt, x) be the polynomial defined as: 

j=0,...,n 

Then by Eq. (3.2), one has p{x\, . . . , x t , —1) = 0. Moreover, for every x > 0, by 
Eq. and (3.1), one has: 

F(xi,...,x t ,x) = Y I a i ' X! XJ = X! a,jPj(x) =p(xi,...,x t ,x). 

j=0,...,n V A=0,...,x / j'=0,...,n 

The proof is thus complete. □ 
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Definition 3 A map F : N — > N is said to be a quasi-polynomial if there 
exist d € N, d > 1, and a family of polynomials in t variables with rational 
coefficients: 

{P(d u d 2 ,-,d t ) I d u ...,dt GN, 0<di <d}, 

where, for every (xi,...,x t ) G N*, if di is the remainder of the division of Xi by 
d, one has: 

F(x l7 ...,x t )= P( dl A 2 ,--.d t ){xi, x t ). 
The number d is called the period of F. 

To simplify the notation, the polynomial P(d 1 ,d 2 ,—,d t ) is denoted Pd x d 2 ---d t - 

Definition 4 Let F : N* — ► N be a map. Given a subset C of N* , F is said 
to be a quasi-polynomial over C if there exists a quasi-polynomial q, such that 
F(x) = q(x), for any x E C . 

Lemma 2 The sum of a finite family of quasi-polynomials is a quasi-polynomial. 

Proof. It suffices to prove the claim for two quasi-polynomials. Let fi,f 2 ■ 
N* — ► N be quasi-polynomials of periods d\ , di respectively and let 

{p ai - a t | Vi = 0, . . . ,t, < a,i < di - 1}, and 

{q bl ... bt | Vi = 0,...,t, 0<bi <d 2 -l} 

be the families of polynomials that define f\ and ,f 2 respectively. Define a new 
quasi-polynomial / as follows. Take d = did 2 as the period of / and, for every 
(ci, . . . , ct) e {0, 1, . . . , d - 1}*, take 

Icf-ct = Va\---at + Qbi-bti 

where, for any i = 1, . . . , t, a{ and bi are the remainders of the division of Cj by 
d\ and d 2 respectively. It is easily checked that the quasi-polynomial / is the 
sum of /i and / 2 . Indeed, if x = (x\, . . . , x t ) G N* and, for every i — 1, . . . , t, 
Xi = a mod d, then one has 

Ci = en mod d\ Xi = di mod d\ 

Ci = bi mod c?2 ^=^> Xi = bi mod d 2 . 
Therefore, if x = (x\ , . . . , x t ) £ N* and Xi = Ci mod d, then we have: 

f(x)= f Cl -c t (x) = P ai -a t (x) + q^-btix) = f\{x) + f 2 (x). 

The claim is thus proved. □ 

Lemma 3 Let F : N — > N be a map, d be a positive integer, and C be a subset 
o/N*. If there exists a family of quasi-polynomials {Fd 1 d 2 -d t I d\,...,d t € 
N, < di < d}, such that , for every (x\, . . . , x t ) € C, with Xi = di mod d, one 
has: F(xi, . . . ,xt) — Fd 1 d 2 ---d t {xi, . . . ,xt), then F is a quasi-polynomial over C . 
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Proof. Let k be the least common multiple of d and of the periods of the quasi- 
polynomials of the set {Fd 1 d 2 -d t | d\, . . . , d t G N, < di < d}. Let (r*i, . . . ,r t ) 
be a tuple of {0, 1, ... k — 1}* and let (xi, . . . ,x t ) G C be such that, for every 
i = 1, t, Xi =ri mod k. Then one can check that F(x±, . . . , Xt) = q(x±, . . . , Xt) 
where q is a polynomial uniquely determined by (ri, . . . , r t ). Indeed, one can 
first observe that the tuple (n, . . . , r t ) uniquely determines, for every i = 1, t, 
the remainder ci^ of the division of Xi by <i since di = r, mode?. By hypothesis, 
one has F(xi, . . . ,Xt) — Fd 1 d 2 ---d t { x ii ■ ■ ■ , ait)- Since k is a multiple of the pe- 
riod of the quasi-polynomial Fd 1 d 2 ---d t , the tuple (r%, . .. ,r t ) also determines a 
polynomial q in the family of polynomials associated with Fd^—dt^ sucn that 
F(xi, ...,x t ) = F dl d 2 ...d t (xi, ...,x t ) = q(xi, ...,x t ). The proof is thus complete. 

□ 

Lemma 4 Let A : N* — > Q be a map such that, for any (xi, . . . , xt) £ N', 

A(xi, . . . , x t ) — bixi H h b t x t , 

where b\, . . . , bt are given rational coefficients. Let C, C be subsets o/N* and let 
a\, . . . , at be non negative integers such that the following properties are satisfied: 
for any [x\, . . . ,Xt) € C, one has 

• X(xi,. ..,x t )>0, 

• t/AsN and X < X(x±, . . . , Xt), then [x\ — Xa±, x-2 — Aa2, . . . , x t — Xa t ) G C" . 
Let p be a quasi-polynomial over C and define the map F as: 

F(x!,...,x t )= ^2 p(xi - Xai,x 2 - Xa 2 , ■ ■ ■ ,x t - Xa t ). 

< A < \(x 1 ,...,x t ) 

Then F is a quasi-polynomial over C . 

Proof. Let d > 1 be the period of the quasi-polynomial p. Let Pdxd 2 ---d t , with 
< di < d — 1, be the polynomials defining p. Consider the set of integers fi: 

0<M < \X(xi,...,x t )] -1, 

and consider on it the partition: 

Fo(xt, ...,x t )U F 1 (xi, . . . ,x t ) U ■ ■ • U Fd-i(xi, . . - ,Xt), (4) 

defined as: for any j = 0, . . . , d — 1: 

H G Fj(xi, . . . , x t ) < /i < \X{x\, . . . , — 1 and \i = j (mod d). 

By Eq. ((!]), for every (ii, . . . , xt) G C, we have: 

d-l 

F(x 1 ,...,x t ) = ^2S j (x 1 ,...,x t ), (5) 

3=0 



where, for any j = 0, . . . , d — 1: 

Sj(xi,...,x t ) = ^ p(xi- fiai,X2- fi02,...,x t - fmt). 

Now we prove that, for any j = 0, . . . , d — 1, Sj is a quasi-polynomial over C. 

Let us fix a tuple (di, . . . , d t ) G {0, 1, . . .d — 1}'. Let {x\, . . . , x t ) be such 
that Xi = di mode?. Then for any [i G Fj(x\, . . . , Xt) one has x\ — fiat = di — 
jdimodd. Now set q — p C ic 2 -c t , where (ci,...,Ct) G {0, 1}' and 

Ci = di — jai mod c?. One has 

Sj(xi,...,x t ) = ^ q{xi~ Liai,x 2 ~ Lia 2 ,.-.,x t - Lia t ). 

fJ,£Fj(x x t ) 



On the other side, by Eq. (|U), one easily checks: 



F j {xi,...,x t )= J + ^eN | 0</i< 



rA(xi, . . .,x t y] - i- j 



It is important to remark that, in the formula above, if X(xi, . . . , Xt) — 0, then 

\X(x u ...,x t )~\ - 1 - j Q 
d 

In this case, since for any < j < d — 1, | — 1 — j\ < d, one has 

-i-i 

Hence A(xi, . . . , xt) = implies that Fj(x±, . . . , x*) is the empty set and thus 
is the null map. Let us consider the map 

S-.N* X {{-1}UN} — >N, 
where 5'(xi, . . . , Xt, x) is defined as: 

X^=o 9(xi - (j + i±d)ai,X2 - (j + fJ.d)a 2 , ...,x t -(j + ^d)a t ) x > 
x = -l. 

By the definition of the map S, for every (xi , . . . , x t ) G C, with x, = mod d, 
one has: 



Sj(xi, . . . ,x t ) = S xi, . . . ,x t 



[A(xi, . . . ,x t )l - 1 - j 
d 



(6) 



Therefore, by applying Lemma[T]to S, there exists a polynomial Q{x\, . . . , xt, x) 
such that, for any (xi, . . . , x t , x) G N* x {{-1} U N}, 



S(xi, . . . ,x t ,x) = Q(xi, . . . ,x t ,x), 



(7) 
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hence, from Eqs. and (0, for every {x\, . . . ,Xt) G C, with Xi = dimodd, 
one has 



Sj{xi, ...,x t ) = Q(xx 



,x t , 



\X(xi, . . .,x t )] -1-3 



)• 



By (c/ 0], Appendix), one has that: 

\\(xi, . . .,x t )] - l-j 



(8) 



(9) 



is a quasi-polynomial in the variables x\, . . . ,Xt- Therefore, by Eq. ([5]) , Sj 
coincides with a quasi-polynomial on the set of points (x\, . . . ,Xt) G C, with 
Xi = di modd. Obviously this fact holds for any (dt, . ■ ■ ,dt) € {0, 1, ... d — 1}* 
and, by Lemma EJ Sj is a quasi-polynomial over C. Finally, the fact that F is 
a quasi-polynomial over C follows from Eq. ([5]) by using Lemma [21 □ 

Lemma 5 Let A : N* — > Q be a map such that, for any (xi, . . . , xt) G N', 

X(xi, ...,x t ) = bixi H h b t x t , 

where b\, . . . , bt are given rational coefficients. 

Let C,C be subsets of N* and let a\,...,at be non negative integers such 
that the following properties are satisfied: for any (x\, . . . ,xt) S C, one has 

• X(xi, ...,x t )>0, 

• for any A G N such that A < X(xi, . . . , x t ), (x\ — Xai, X2 — Aa2, . . . , x t — 
Xa t ) 6 C. 

Let p be a quasi-polynomial over C and define the map F as: 

F(xi,...,x t ) = ^2 p{xi - Xai,x 2 - Xa 2 , ■ ■ ■ ,x t - Xa t ). 

< A < \(xi,...,x t ) 



Then F is a quasi-polynomial over C . 

Proof. The proof of Lemma [5] is the same of that of Lemma 2] except the point 
we describe now. In the sum above that defines the map F, the index A runs 
over the set of integers of the closed interval [0, X(xi, . . . , Xt)] so that: 



A < [X(x 1 



,x t )\. 



Therefore, in order to prove the claim, one has to prove a slightly modified 
version of Eq. ([9]) of LemmalU that is: for any j = 0, . . . , d — 1, 

|_A(a;i, . . .,x t )\ - j 



is a quasi-polynomial with rational coefficients in the variables x\, . . . ,x t . This 
can be done by using an argument very similar to that one adopted in the proof 
ofEq. ©. □ 
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Lemma 6 Let Xi, X2 '■ N* — ► Q be 2 maps such that, for any (xi, . . . , xt) G N', 

Ai(xi, ...,x t ) = bixi H h b t x t , X 2 {x 1 , ...,x t ) = c x x x ~\ h c t x t 

where b\ , . . . , bt and c\ , . . . , Ct are given rational coefficients. 

Let C be a subset ofW and let 01, . . . , at be non negative integers. Suppose 
that, for any (xi, . . . , xt) S C , one has: 

< Ai(xi, . . . ,x t ) < A 2 (xi, . . .,x t ), 

and, for any A G N such that Ai(xi, . . . ,Xt) < A < X2(xi, . . . ,Xt), one has: 

(xi - Aai, x 2 - Aa 2 , ...,x t - Xa t ) e C", 

where C is a given subset o/N . Let p be a quasi-polynomial over C and define 
the map F as: 

F(xi,...,x t ) = ^ p(xi - Aai,x 2 - Aa 2 , .. .,x t - Aa 4 ), 

Ai < A < A 2 

where Ai = Ai(xi, . . . , x t ) and A 2 = A 2 (xi, . . . ,xt). Then F is a quasi-polynomial 
over C . The same result holds whenever the index X runs in the set of integers 
of the intervals: 

(Ai,A 2 ), (Ai,A 2 ], [Ai,A 2 ). 
Proof. Let us solve the case when A runs in the interval [Ai, A 2 ]. Write 

F(xx, ...,x t ) = S 1 (x 1 , . . . ,x t ) - S 2 (x 1 , . . .,x t ), (10) 

where: 

Si(xi,...,x t ) = ^2 p(xi - Aai,x 2 - Aa 2 , . . . ,x t - Xa t ), 
< a < A 2 

and 

5 2 (xi, . . . , x t ) = p{x\ - Aai,x 2 - Aa 2 , . . . ,x t - Xa t ). 

0<A<Ai 

By applying Lemma O to S\ and Lemma 0] to S* 2 , we have that S\ and 5 2 are 
quasi-polynomial and by Lemma [2l so is S% — 5 2 . The claim now follows from 
Eq. (fT0|) . The other three cases are similarly proved. □ 

Lemma 7 Assuming the same hypotheses of Lemma\^ the function 

S(xi, . . . ,x t ) = ^2 p(xi - Aai,x 2 - Aa 2 , . . . ,x t - Xa t ). 

X(x±,...,xt) < A < X(x± , . . .,xt ) 

is a quasi-polynomial over C. 
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Proof. It is a direct consequence of Lemma HI assuming Xi(x\, . . . , Xt) = 
X 2 (xi,...,x t ) = X(xi,...,x t ) □ 



Now we want to define some suitable regions of R*. More precisely, our 
regions will be polyhedral cones determined by a family of hyperplanes passing 
through the origin. We proceed as follows. Let it be a plane of R*. Let us fix 
an equation for 7r denoted by 7r(x) = 0. We associate with n a map 

U :R* — > {+,-,0} 

defined as: for any x G R', 

!+ if tt{x) > 0, 
if 7r(^) = 0, 
- ifvr(a;)<0. 

We remark that the map defined above depends upon the plane ir and its equa- 
tion in the obvious geometrical way. We can now give the following important 
two definitions. 

Definition 5 Let H = {tt\, . . . , 7r m } be a family of planes o/R* that satisfy the 
following property: 

• II includes the coordinate planes, that is, the planes defined by the equa- 
tions X(. = 0, £ — 0, . . . , t; 

• every plane of II passes through the origin. 

Let ~ be the equivalence defined over the set IR' as: for any x, x' G R , 

x ~ x' t=> Vi = l,...,m, Ui(x) = U t (x'). 
A subset C of R* is called a region (with respect to II) if it is a coset of ~ . 

It may be useful to keep in mind that the singleton composed by the origin 
is a region. Moreover if t = 2, the set of all points of R* \ {0} of every line of LT 
is a region also. 

Definition 6 Let F : N* — ► N be a map. Then F is said to be a box spline 
in N* if there exists a partition C — {C±, . . . , C y } of regions of N* - defined by 
a family of planes satisfying Definition [3| - and a family p±, . . . ,p y of quasi- 
polynomials, every one of which is associated with exactly a region ofC, such 
that, for any (xi, . . . , a;*) G N', one has: 

F(xi, ...,x t ) =p a (xi, . . -,x t ), 

where a is the index of the region C a that contains (xi, . . . , Xt). 

Lemma 8 The sum of a finite family of box splines is a box spline. 
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Proof. It suffices to prove the claim for two box splines. Let F\ and F2 be two 
box splines and let C — {C\, . . . , C y } and V = {Z?i, . . . , D z } be the families of 
regions of Fi and F2 respectively. Moreover, let {pi, . . . ,p y } and {qi, . . . ,q z } 
be the families of quasi polynomials of F\ and F2 respectively. 

Consider the box spline defined as follows. Let £ be the partition of regions 
of N* given by the intersection of C and V respectively. It is worth noticing that 
£ is determined by the union of the two families of hyperplanes that define C 
and V respectively. Then we associate the map r; m = pi+ q m with every region 
Ei m of £ . By Lemma [21 r; m is a quasi-polynomial. For any x G N* we have 

Fi(x) =pi(x), F 2 (x) = q m (x), 

where / and m are the indices of the regions Ci and D m that contain x. Hence 
we have 

Fi(x)+F 2 (x) = r lm (x), 

while x belongs to the region Ei m . Since r; m is the quasi polynomial associated 
with Ei m , this proves that F\ + F2 is equal to the box spline defined above. □ 

Lemma 9 Let F : N* 1+ * 2 — > N be a box spline in N* 1+ ' 2 . Then there ex- 
ists a map G : N* 1 — > N, which is a box spline in N 1 , such that, for every 
(xx, . . . , Xt t ) G N* 1 , the following equality holds: 
G(x\, . . . ,x tl ) = F(xi, . . . ,x tl , hxu hxi ). 

ti-times 

Proof. Let Uf = {fti, ■ ■ ■ , 7r m } be the family of planes of R tl+t2 associated with 
F. 

Recall that Hp satisfies the following properties: 

• Ili? includes the coordinate planes, that is, the planes defined by the equa- 
tions xe — 0, I = 0, . . . , t\ + t<i\ 

• every plane of IIf passes through the origin. 

Now let 7r(:Ei, . . . , xt 1 +t 2 ) = J2i=i t ± +t 2 @ iXi — be a plane in 11^. Then the 
plane: 

■K'{xi,...,X tl ) =Et=l,...,t I Pi^i +J2j=t 1 + l,...,t 1+ t 2 Pjhxx=0 

is a plane of M* 1 through the origin. We define as IIg the family of all such 
planes. It is obvious that all coordinates planes belong to lie. 

Let now x, — (x%, ...,5^) be a point in N* 1 . Since F is a box-spline, F 
associates to the point x* — (xi, x^, hxi, hx\) a unique quasi-polynomial 

^2-times 

p(x). 

We want to show that the region of x w.r.t. the planes in ITg determines 
p(x) univocally. Let it = E»=i t +t @ iXi = be a plane in lip. Let ir' be the 
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corresponding plane of Hq- Then it is obvious that tt(x*) > (or ir(x*) = 0, or 
7r(x*) < 0) if and only if ir'(x) > (or, respectively, ir'(x) = 0, or, respectively, 
tt'(x) < 0). Then the position of x w.r.t. He determines p(x) univocally. 

Now, let d the period of p{x). Observe that a given point x = (xi, Xt 1 ) in 
N* 1 , gives rise to the remainders (d\, . . . , d tl , d\, dl) modulo d, where: 

i 2 -times 

dl = hxi mod(rf). 

It follows that, each sequence of remainders (d\, . . . , d tl ) modulo d specifies 
the unique polynomial P{di,di,—,dt ,d{ d\) °f the quasi-polynomial p(x). 

Therefore we let correspond to p(x) the quasi-polynomial q(x) : N* 1 — ► N, 
which has period d and to each sequence of of remainders (di, . . . , d tl ) modulo 
d associates the polynomial q(d 1 ,d 2 ,---.d tl ) '■ N* 1 — > with rational coefficients, 
defined by; 



Q , (d 1 ,d 2 ,-,dt 1 )( a;i ' —■> x ti)= P(d- L ,d^,-,dt 1 ,dl,-,dl){ x ^ x *i ' hx i) 

where (xi, x^) is such that di = Xi mod(d) , for every i — 1, ti, and 
d\ = hx\ mod(d). 

The box-spline G is therefore completely specified and this ends the proof. 

□ 

The following lemma is a crucial tool in the proof of the main result of this 
section. 

Lemma 10 Let G : N* — ► N be a box spline and let a\,...,a t E N with 
(a\,...,a t ) ^ (0, ...,0). Consider the map A : N* — ► N that associates with 
every (xi, . . . , Xt) G N* ; the value 

A(xi,...,x t ) =min| ^ | a» ^ o| . 

Let S : N* — ► N be the map defined as: for every (x\, . . . , x t ) € N*, 

S(xi,...,x t )= ^ G(xi - \ai,x 2 - \a 2 ,...,x t - \a t ). (11) 

0<A<A(xi,...,x t ) 

Then S is a box spline. 

Proof. In order to prove the claim, we first associate with the map S a new 
family of regions that we define now. Let II be the family of planes associated 
with the box spline G. For any (x\, x 2 , ■ ■ ■ , x t ) £ N* consider the line defined 
by the equation parameterized by A: 

(ari - Xai,x 2 - Aa 2 , . . . ,x t - Aa t ). (12) 
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Let 7r be a plane of the family II and let tt(x) — J2i=i t @iXi = be its 
equation. The value of A that defines the point of meeting of the line (| 1 2[) with 
7r is easily computed. Indeed, A is such that 



0i(xi - Aaj) = 0, 



«=i,...,t 



so that 



(13) 



i=l t t=l t 



which gives 



A 



i=l,..., 



E 




(14) 



where 



7 = 2J Act- 



»=i,...,t 

It is worth to remark that Eq. (|14p is not defined whenever 



7 = E = 0. 



(15) 



»=i,...,t 



Let us first treat Eq. (|15[1 . Here, either the line of Eq. (fl2|) belongs to 7r or such 
a line is parallel to n. Therefore, for every point x of the line of Eq. (fT2|) . the 
value of /tt(x) is constant so that 7r is not relevant in determining a change of 
region when a point is moving on the line of Eq. (|12p . Because of this remark, 
we shall consider only planes of II for which Eq. (fl~5|) does not hold. Denote IT' 
this set of planes. For any 7r £ II', with equation -k(x) = 2j=i t Pi x i = 0; 
consider the homogeneous linear polynomial 



where 7 = J2i=i t @i a i- We remark that for any (xi, X2, ■ ■ ■ , %t) € N*, the line 
parameterized by Eq. (|12|) meets the plane 7r in the point corresponding to the 
parameter X — X 7T (xi, . . . ,x t ). 

Consider an enumeration of the planes of the set II and denote by < the 
linear order on II defined by such enumeration. Consider the new family II of 
planes defined by the following sets of equations: 

1. n(x) = 0, 7T e n 

2. Xm!-'(xi, . . . ,Xk) = 0, with 7r,7r' e II', tt < n', and A,™-' {xi , ■ . ■ , Xk) = 
X n (xi,..., x k ) - A ff '(xi, . . . , Xk). 

Call C the family of regions of N* defined by II. 
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We now associate with every region of C a quasi-polynomial. In order to do this, 
we need to establish some preliminary facts. Let us fix now a region C of C and 
let x = (x\, . . . , Xt) be a point of N* that belongs to C. Let i be such that 



A <*>=S- 



Observe that, for any other point x' = (x[ , . . . , x' t ) in C, one has 

A(x') - ^. 

(Xi 



Indeed, it is enough to prove that, for any given pair of distinct indices 
we have: 

— < f2. x j < X J 

This is equivalent to say that: 

A 7rTr /(xi, . . . ,x t ) < A^/(xi, . . . ,x' t ) < 0, 

where it, tt' are the planes Xi = and Xj = respectively. The previous equiva- 
lence is true because x and x' belong to the same region of C. 

Another important fact is the following. Let us consider any point x of the 
region C of C. Consider the subset of planes of II': 

{tti, . . . , 7r m } = {tt g n' I < A^Cx) < A(x)}. 

We can always assume, possibly changing the enumeration of the above planes, 
that 

< (x) < • • • < K m (x) < A(x). 
Remark. Observe that, for any other point x' of C, one has 

{tti, . . . , 7r m } = {tt e n' | < A^(x') < A(x)}. 

and 

< Xm(x') <■■■< K m {x') < A(x'). 

The remark above can be proved by using an argument very similar to that 
used to prove the previous condition. We suppose that the above inequalities 
are strict, i.e. < A 7ri (x) < • • • < A^^x) < A(x). In this before, one 

proves that the same inequalities are strict for any other point x' of the region 
C. The case when the inequalities are not all strict can be treated similarly. 

From now on, by the sake of clarity, for any x = (xi,...,x t ) s N we set 
y\(x) = (xi - Aai, ...,x t - Xat). 
Consider the following sets: 
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. Y (x) = {y x (x) | AeNn[0,V(x))}, 

• Y m (x) = {y x (x) | AeNn(A Wm ,A(x))}, 

• Yi{x) = {y x (x) | AeNn (A^ (x),X ni+1 (x))}, i = 1, . . . ,m — 1. 

• Zi{x) = {yxix) | AeNn{A Ij (x)}},j = l 1 ...,m. 

• Z m+ i(ar) = {y A (a;) | A eNn{A(i)}}. 

We are now able to associate a quasi-polynomial with the region C of C. For 
this purpose, take two points x,x' in C. By the facts discussed before, one 
has that the lines of Eq. (p~2|) associated with x and a;' respectively, meet the 
planes of II' in the same order. We recall, that a change of region on the generic 
line of Eq. (TT2]) happens only when the line meets a plane of II'. Therefore, 
since C is a refinement of C and x and x' are in a same region with respect 
to C, the above conditions imply that, for every i — 0, . . . , m, the two sets of 
points Yi(x) and Yi(x') are subsets of a same common region of C. Hence there 
exists a quasi-polynomial pi, depending on i and on the region C, such that, 
for any y E Yi(x) and for any y' E Fj(x'), G(y) = pi(y), G{y') = Pi(y'). By the 
previous remark and by Lemma[51 one has that, for any i = 0, m, there exists 
a quasi-polynomial qi, depending on i, and on C, such that for any x 6 C 

qi ( x )= J2 °(y)- 

yeYi(x) 

Observe that, since x and x' are in the same region C, as before one derives 
that Zi(x) and Zi(x') are in the same region with respect to C. Therefore, as 
before, by applying Lemma [7] there exists a quasi-polynomial r^, depending on 
i and on C, such that for any x E C 

n{x)= G (y)- 

yez,(x) 

On the other hand, by Eq. (jTTj) . we have that, for any (x%, . . . , xt) E C, 
S(xi, . . . , Xt) is equal to: 

q {x) +ri(x) +qi(x) + r 2 (x) + q 2 {x) H r m (x) + q m {x) +r m+1 (x). (16) 

Thus, iS(a;i , . . . , Xt) on the region C is represented as a sum of quasi-polynomials. 
This, together with Lemma applied to Eq. (|16p imply that the map S is a 
quasi-polynomial over every region of C. The proof of the claim is thus complete. 

□ 

Theorem 2 Let 

S : N* — ► N 
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be the function which counts, for any vector (n\, . . . ,n t ) 6 N*, the number of 
distinct non-negative solutions of a given Diophantine system: 

cinxi + a 12 x 2 + ■ ■ ■ + aikXk = n\ 
a 2 \Xi + a 22 x 2 H h a 2k x k = n 2 

[ [ (17) 

a t \Xx + a t2 x 2 H h citkXk = 

where the numbers S N and, for every i = 1, . . . , k, there exists j = 1, . . . , t 
such that aij ^ 0. The function S is a box spline. Moreover such box spline can 
be effectively constructed starting from the coefficients of the system. 

Proof. For any vector (m, . . . , ni) 6 N* , let Sol(ni, . . . , n t ) be the set of the non- 
negative solutions of the Diophantine system (JTTJ) and denote by S : N* — > N, 
the map defined as: for any vector (m, . . . , n t ) G N*, 

S{m, ...,n t ) = Card(Sol(ni, . . . ,n t )), 

that is, it associates with every vector (ni, . . . , n t ) the number of non-negative 
distinct solutions of the system (jTTj) . Let us prove that the map S is a box 
spline. For this purpose, we proceed by induction on the number of unknowns 
of the system (JTTJ). We start by proving the basis of the induction. In this case, 
our system has one unknown, say x, and it can be written as: 

a\x = ni 
a 2 x = n 2 

< 

a t x = n t 

The system has solutions (and, in this case, it is unique) if and only if there 
exists A 6 N such that: 

X(a%, ...,a t ) = (Aoi, . . . , Xa t ) = (m, . . . ,n t ). (18) 

Let us consider the line £ (through the origin) defined by the parametric equation 
(fT5|) . The line £ can be determined as the intersection of suitable planes through 
the origin. Let us consider the family of regions defined by the set of these planes 
together with the coordinate planes. One can easily associate with every region 
a quasipolynomial. For this purpose, we remark that the set of points of the 
line I with integral coordinates, without the origin, is a region. On this region, 
the counting function of the system takes the value or 1. Therefore this map 
coincides with the quasi-polynomial given by p — 0, q = 1 with the periodical 
rule d = lcm{ai 1 . . . , at}. To any other region, we associate p. The basis of the 
induction is thus proved. 
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Let us now prove the inductive step. If (x\, . . . , Xk) £ Sol(ni, . . . , n t ), the 
system (fT7|) can be written as: 



012X2 

0-22X2 



a t 2X2 



aikXk 
a.2kXk 



atkXk 



rii — ai\X\ 
n 2 - a 2 ixi 



n t - at\X\. 



(19) 



This implies that: 

m - CLuXi > 0, ri2 — a 2 \xi > 0, n t - a n xi > 0, 
so that, since x\ must be an integer > 0, one has: 



"l 



"2 



< xi < — , < xi < 

an o 2 i 



< xi < 



an ' 



and thus: 



< xi < A(ni, . . . ,n t ), 
where the map A : N* — > N is defined as: 



A(xi, . . . ,x t ) = mini — | an ^ 

an 



We remark that, since the vector (an, asi, . . . , aji) ^ (0, 0, . . 
A is well defined. Set K — [A(xi, . . . ,xt)\. We can write Sol(ni, 



(20) 

0), the map 
. . , n t ) as: 



Sol(ni,...,n t ) = (0 X Sol ) U (1 X Soli) U ... U (KxSol K ), 



(21) 



where, for every i = 0, . . . , K, Sol^ denotes the set of non-negative solutions of 
the Diophantine system: 



Ol2X2 
022X 2 



aikXk 
a2kXk 



n\ — am 
n 2 - a 2 ii 



(22) 



a t2 x 2 H h a tk x k = n t - a tx i. 

By Eq. (I2ip . for any (n 1: . . . , n t ) e N*, we have: 

5(ni,...,n t ) = Card(Soh). 



(23) 



i=0,...,if 



By applying the inductive hypothesis to the system (12"2")) , we have that there 
exists a box spline G : N* — > N such that, for any [n\, . . . , n t ) € N', if < i < 
K, 

Card(Soli) = G(nx - aui, n 2 - a 2 ii, . . . ,n t - atxi), (24) 
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so that, by Eq. (23]) and Eq. (J2JJ, one has: 



S(nx,...,nt)= ^ G(ni - Aan,n 2 - Aa 2 i, 

0<A<A(a:i,...,an t ) 



,rt f -Aa t i). (25) 



By Eq. (|2"5|) , the fact that S is a box spline follows from Lemma ITD1 Finally we 
remark that the proof gives an effective procedure to construct the claimed box 
spline that describes the map <S. □ 



Corollary 1 Let 



S : N* 



N 



be the function which counts, for any vector (n\, . . . ,n t ) £ N* ; the number of 
distinct non-negative solutions of a given Diophantine system: 



aio + an^i + «12^2 + 

O20 + 0-21X1 + a22X2 + 



+ aikx k = m 
+ a2kXk = n 2 



(26) 



am + a t ixi + a t2 x 2 H V a t kXk = n t , 



where the numbers atj S N and, for every i = 1, . .., k, there exists j = l,...,t 
such that aij ^ 0. Set ao = (aio, • • ■ , ato)- Then, on the set of all vectors 
x £ N* with x > ao, the map S is a box spline. Moreover such box spline can be 
effectively constructed starting from the coefficients of the system. 



Proof. First consider the system 



anxi + ai 2 x 2 
a 2 ixi + a 22 x 2 



+ aikXk = ni 
+ a 2 kXk = n 2 



(27) 



a t ixi + a t2 x 2 H h a tk Xk 



According to Theorem [2l there exists a box spline F that counts, for every 
(rii, . . . , nt) S N*, the number of the solutions of the diophantine system ([27]). 
Let C — {Ci, . . . , C s } be the partition of N* in regions and let {pi, . . . ,p s } be 
the family of quasi-polynomials that define F. 

Let ao = (aio, ■ ■ ■ , (Ho) be the vector whose components are the entries of 
the first column of the matrix of the system ([2"6"| . For every rj £ N* with r\ > ao, 
the components of the vector r\ — ao are non negative integers so that: 

S(V) =Pj(v - a ), 

where j is the index of the region of the family C that contains the vector r\ — oq. 
This concludes the proof. □ 
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Example 1 For the sake of clarity, we find useful to show the proof of Theorem 
[2] on the following very simple example. Consider the Diophantine system: 

f xi + 2x 2 = ni , 28 -. 
\ 2xi + 3x 2 = n 2 , 

where n\, ri2 £ N and let F : N 2 — > N be the counting function of the system 
(|2"8|) . By following the proof of Theorem^ together with that of Lemma [TU1 we 
construct the partition of N 2 in regions and the family of quasi-polynomials that 
describe the function F. In the sequel, the following notation is adopted: [x\ , x%) 
and (n 1,712) are respectively the vector of the unknowns and the vector of the 
non homogeneous terms of the system, while x, y are free variables over the set 
N. Observe that x\,X2 gives a solution of ([28]) if and only if (ni — x±, 712 — 2xi) — 
(2i,3t), t > 0. 

Therefore consider the Diophantine system: 

( l X 2 = ni (29) 
[ 3x 2 = n 2 

where ni,7i2 £ N. Let G : N 2 — > N be the counting function of the system 
(|2"5|) . Let LI = {tti, 7T2, ^3} be the set of the lines defined by the equations: 

ni(x,y)=x = 0, n 2 (x,y) = y = 0, ir 3 (x, y) = 3x - 2y = 0. 

Let 1Z be the partition of N 2 determined by II. Then 1Z is formed by the fol- 
lowing 6 regions: 

R = {(0,0)}, i?i = {(x,Q) : x > 0}, R 2 = {(x,y) : 3x > 2y, x,y > 0}, 
#3 = {{z,y) ■ 3x = 2y, x,y > 0}, R 4 = {(x,y) : 3x < 2y, x, y > 0}, 
R 5 = {(0,y) : y>0}. 

Let V be the family of polynomials given by: 

Po(x,v) =P3(x,y) = 1, Pi(x,y) =p 2 (x,y) =p 4 (x,y) =p 5 (x,y) = 0. 

One can check that the box spline determined by 1Z and V is the function 
G. 

Let (x, y) be a given point of N 2 and let £(X) be the line represented by the 
equation parameterized by A: 

(x - A, y - A). 

For every i = 1,2,3, let X ni (x,y) be the value of A that defines the point of 
meeting of the line £(X) with 7Tj. Then one has: 

V (x, y) = x, X W2 (x, y) = y/2, A W3 (x, y)=2y- 3x. 

Consider the new family LI of lines defined by the following sets of equations: 
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1. TTi(x, y) = 0, i = 1, 2, 3, (that is, the lines in II), 

2. For every pair i,j of indices with 1 < i < j < 3, let 

{x, y) = Ki (x, y) - K 3 {x, y) = 0. 

It easily checked that there exists exactly one line of the previous form (2) which 
is represented by the equation y — 2x = 0. Thus one has: 

fi = {x = 0, y = 0, 2x - 3y = 0, 2x - y = 0}. 
If we denote by 7?. the family of regions of N 2 defined by II, we have: 
R = {(0,0)}, i?i = {(x,0) : x > 0}, R 2 = {(x,y) : 3x > 2y, y > 0}, 
#3 = {(x, y) : ix = 2y, y > 0}, R 4 = {(x, y) : 3x < 2y, 2x > y, x,y > 0}, 
R5 = {(x,y) : 2x = y, y> 0}, i? 6 = {(x, y) : 2x < y, x, y > 0}, 
R 7 = {(0,y) : y>0}. 

Now we associate with every region Ri of R a quasi-polynomial q{. Actually, 
set: 

Qo = 93 = 94 = 95 = 1 and & = g 2 = % = <?7 = 0. 

One can check that the box spline determined by 1Z together with the list of 
polynomials above is the function F. 

4 The growth function of a semi-linear set 

The main result of this section is that the growth function of a semi-linear set 
is a box-spline. Let us start by introducing a definition. 

Definition 7 Let X be a subset of Z* and let t\,t 2 G N such that t = t\ +t 2 - 
We associate with X a function 

Q + x, tl>ta :N* — >N 
that returns, for every (n\, . . . ,n tl ,mi . . . ,m t2 ) £ N* 7 the number of elements 
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(x±, . . . , Xt) € X such that: 



X\ 



(30) 



< Xt-t+i < mi 
< x t < m t , . 



The function Q + x,ti,t 2 * s called the (restricted) generalized growth function of 
X. 

Theorem 3 Let X be a semi-linear set o/N* and let t\,t2 E N such that t — 
t\ +t%- Then G + x,t u t 2 * s a oox spline. 

Proof. Let ti,ts S N be given as above and let X be a semi-linear set of N*. To 
avoid a heavy notation, in this proof, we suppress the dependency of G + x on 
ti and t2 and we simply write Qx- By Theorem [TJ X is semi-simple so that 



X= |J Xi 



i=l,. 



is a finite and disjoint union of simple sets of N*. As a straightforward conse- 
quence, one can easily check that, for every (n%, n t ) S N*, 



Qx(ni,—,n t )= ^ Gxi(ni,—, 



m ) 



In order to obtain the claim, by the equality above and by Lemma [8] it is enough 
to prove that, for every simple set Y of N*, Qy is a box spline. 
For this purpose, let 

Y = b + bf + --- + bf, 
be a representation of Y as a simple set where, for every j = 0, k, 

bj = {b jl ,---,b jt )€N t . 
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Therefore we can write Eq. ([50)) for the set Y as a system of t inequalities in k 
unkowns yi,...,y k : 



boi 
b 2 



buyi 



buyi 

^222/2 



bkiyk 
bk2Vk 



boa + biaVi + b 2tl V2 H V b ktl y k 

bott+i + hti+iVi + b 2 t l +iV2 + ■ 
&o*i+2 + b Ul+2 yi + b 2tl +2y2 + ■ 



hi 

n 2 



htVi + b2ty2 



bktVk 



b ktl+ iy k < mi 
bkt!+2yk < m 2 



< m t2 . 



(31) 



Since Y is a simple set of N*, there exists a bijection between the set of non- 
negative solutions of the previous system and the set of elements of Y that 
satisfy Eq. ([50)1 . Consider now the Diophantine system of equations obtained 
from that of Eq. (|3ip where z\, z 2 , ■ z t2 form a set of t unknowns disjoint from 
the set of unknowns yi,...,yk- 



boi + b n yi + b 2 \y2 H V b kl y k = m 

bo2 + &12J/1 + 6222/2 H h b k 2yk = n 2 

bo tl + b ltl yi + b 2tl y 2 -\ \-b ktl y k =n tl 

(32) 

bou+i + ha+iVi + b 2tl +iy2 H 1- b ktl+1 y k + z\ =m 1 

froti+2 + b ltl+2 yi + b 2tl +2y2 H h b ktl+2 y k + z 2 = m 2 

&0* + feitj/i + &2t2/2 H ^b kt y k + z t2 =m t2 . 



One can easily check that the number of non-negative solutions of the Dio- 
phantine system of inequalities of Eq. (|3ip is equal to the number of non- 
negative solutions of the Diophantine system of Eq. (f3"2")l . By Corollary [1] 
applied to the latter system, we have that its the counting function is a box 
spline. By the previous facts, the function Qy is a box spline as well and this 
concludes the proof. □ 

Now we consider a first extension of Theorem [3] to semi-linear sets over Z*. 
For this purpose the following lemma is useful (see [2]). 
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Lemma 11 If X is linear in Z* then InN' is semi-linear in N*. 

Corollary 2 If X is semi-linear in 1} then XflN' is semi-linear in N . 

Proof. We can write X as a finite union of linear sets X\, Xt oil} . Thus we 
have 

inN'= |J x't, 

i=l,...,i 

where, for every i = 1, ...,£, X[ — X{ flN*. The claim now follows from the fact 
that, by Lemma fill for every i = 1, I, Xi is semi-linear in N*. □ 

Corollary 3 Lei I fc u semi-linear set of Z and /ei ij.,t2 € N swc/i i/ia£ 
£ = t\ +t%. Set X' = X nl . TTien ti,t 2 * s a ^° 2 ' spline. 

Proof. Since, by hypotheses, X is semi- linear in Z*, by Corollary [21 X' is semi- 
linear in N*. The claim follows by applying Theorem [3] to X'. □ 

Definition 8 Let X be a subset of 1} and let ii,i2 G N smc/j f/iat i = ti +i2. 
We associate with X a function 

Qx, t lM ■ N* — ► N 

</ia< returns, for every (n\, . . . ,n tl ,m\ . . . ,m t2 ) G N , i/ie number of elements 
(xi, . . . , Xt) G X such that: 

\x±\ = n\ 



F*il = n *i 



|x t | < m t2 . 

The function Qx, ti t 2 * s called the generalized growth function of X. 

Remark 5 If ti = 0, we obtain the function defined in Eq. (TTJ) of Definition [2] 
that we have called the growth function of X . 

Theorem 4 Let X be a semi-linear set of 1} and let tx, t% G N such that t — 
t\ +t2. Then Qx. t t .t 2 * s a box spline. 
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In order to prove the theorem, we need some preliminaries. First we refresh 
some notions given in Section [2] Let IT be the family of planes 7Tj of equation 
xi = 0, for every i = 1, ...,t and let ~ be the equivalence relation introduced 
in Definition \5\ We recall that a region is a coset of Z* with respect to ~ . 
Let C be the family of regions associated with II. It is easily checked that every 
region C is a semilinear set of Z* . On the other hand, we recall that the family 
of semi-linear sets of Z* is closed with respect to the Boolean set operations. 
These facts allows one to obtain the following useful result. 

Lemma 12 Let X be a semi-linear set ofl} and let C € C. Then the setXdC 
is still semi-linear in Z . 

Lemma 13 Let X be a semi-linear set of Z and let t\,ti £ N such that t = 
t\ +t%. Let C € C and set Y = X D C. There exists a bijection 

f : Z* — > Z\ 

that satisfies the following conditions: 

1. f(Y)CN f , 

2. the map f preserves semi-linearity in 1} , 

3. For every n S N*, Qy, t u tM = G + f{Y), t u tM- 
Therefore, the map Qy, ti,t 2 w a box spline. 

Proof. For the sake of simplicity, we study the case t = 2, the general case being 
treated similarly. Set N+ = N \ {0}. The family C is given by: 

• C x = {(0,0)}, 

• C 2 =N 2 + , 

• C 3 = {(-x,y)}\x, y eE+}, 

• C 4 = -y)} KyeN + }, 
. C 5 = {( Xl -y)}\x,yeE + }, 

• C 6 = {(0,y)}\yeN + }, 

• C 7 = {(0,-y)}|yeN + }, 

. C 8 = {(z,0)}|zeN+}, 
. C 9 = {K0)}| 3 ;eN + }. 
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Set Y — X (~l C, where C G C and X is a semi-linear set of Z 2 . If C = Ci 
then the result is trivial. If C G {C2, C6, Cs}, then the result follows by taking 
/ as the identity on I? and by applying Corollary [3] to G + y, ti,* 2 ( r 7)- Assume 
C = C3 and define the function / : Z 2 — ► Z 2 as: for every (z\, 22) G Z 2 

f{z 1 ,Z2) = (-z 1 ,z 2 ). 

One immediately has /(Y) C N 2 which gives Condition (1). Moreover, the 
function / is an isometry which gives Condition (3). Let us prove Condition 
(2). If X is a linear set of Z', t > 2, and B is a representation of X, then one 
can easily check that f(B) is a representation of f(X). Then Condition (2) 
follows from the fact that a semilinear set is a finite union of linear sets. Finally 
the fact that the function Qy, t t ,t 2 is a box spline follows from Condition (3) 
by applying Corollary [3] to Q + fry), t 1 ,t 2 - We remark that the other two cases 
C4, C5, C7, C9 can be treated similarly as C3. □ 

We are now in position to prove Theorem [H 

Proof of Theorem [4] Let X be a semi- linear set of Z* and let ii, G N 

such that t = £1 + £2 • The set X can be written as a finite and disjoint union 

x= |J (inc). 

cec 

Therefore, for any 77 G N', one has: 

Ox, txtein) = G(*nc, t!,t 2 )0/)- (34) 
cec 

If C is any region of Z', the fact that Gxnc, t ± ,t 2 is a box spline follows from 
Lemma [131 Finally the claim follows by applying Lemma [8] to Eq. (|34|) . □ 

5 A combinatorial proof of a theorem of Dah- 
men and Micchelli. 

Let A be a matrix in Z* xn , with n > t. Assume that the following condition 
holds: 

V X G Z™, X > 0™, AX = 0* X = 0™. (35) 
The following property holds. 

Lemma 14 Lei A k a matrix in Z tx ™ 7 ««t/i n > t 7 satisfying H35\) . If b is a 
vector in Z* ; i/ien the number of non-negative integer solutions of the system 
AX = b is always finite. 
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Proof. By contradiction, assume that the number of non-negative integer solu- 
tions of the system AX = b is infinite. Since N" is well quasi-ordered, there exist 
two solutions X\ and X2 of the system AX = b such that Xi > X2 . Thus the 
vector X\ — X2 is a non null non- negative solution of the system AX = 0* . This 
contradicts the hypothesis that A satisfies (I55|) . The claim is thus proved. □ 

Therefore, given a matrix A satisfying (|35|) . we can define a function 

C A : Z* — > N, 

which associates, with every vector b G Z*, the number of non-negative integer 
solutions of the system AX = b. The following problem has been addressed and 
solved in [5]. 

Theorem 5 (Dahmen and Micchelli,1988) Let A be a matrix in Z* x ", with 
n >t, that satisfies i35\) . Then the function Ca is a box spline in Z*. 

In the following, we give a new proof of the latter theorem. Our proof is based 
on the results of the previous sections and on the properties of semi-simple sets 
in Z*. 

Let SA,b be the set of all non-negative real solutions of the system AX = b. 
with b G Z* . Now, we want to show that Condition ([35]) on the matrix A implies 
that SA,b is a bounded set in R™ . First, as done in [5], we observe that (|3"5|) is 
equivalent to say that the convex hull Ha of the columns et^'s of the matrix A, 
that is, 

Ha = \ (J2 a ^) J= i,...,t 

I i=l 

does not contain the origin 0*. 

Let X be a closed convex set of R* and let 1 be a point of the boundary of X. 
An hyperplane ir such that x G it and X is contained in one of the two closed 
halfspaces determined by n, is called support hyperplane of X at point x. The 
following well-known result of Convex geometry holds. 

Theorem 6 (fltf. Theorem 4-1) Given a closed convex subset X o/R* and a 
point x of the boundary of X, there exists a support hyperplane of X at x, not 
necessarily unique. 

Denote by ||ir|| the Euclidean norm of a point x € R*. Let A be a matrix in 
Z tx ™, with n > t, that satisfies (|35]l. Hence we have 0* ^ Ha- Since Ha is 
a convex and closed set of R', by applying Theorem [5] to Ha, there exists a 
hyperplane tt that separates Ha from the origin 0*. If 5 is the minimal distance 
of 0* from n, then 5 > so that we have: 

Corollary 4 Let A be a matrix in Z* x ™, with n > t, that satisfies &35\) . Then 
there exists a real number S > such that, for every x G Ha, \\x\\ > S. 



Cti G 



i=l 



= 1 
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Proposition 1 Let A be a matrix in Z tx ™ , with n >t, that satisfies (SEj) . If b 
is a vector in 1} , then the set SA,b is bounded in K™ . 

Proof. We argue by contradiction. Assume that is not bounded. Then 
there exists a partition of {1, . . . , n} in two sets / and J such that the following 
property holds: for every i £ / (resp. j G J), the set of values that appear 
in every vector of SA,b at position i (resp. j) is unbounded (resp. bounded). 
Let io be in /. We can define an infinite sequence of solutions (X r ) r >i which 
is unbounded and such that for every r, < X£ + . Consider now any other 
coordinate i\ ^ io in I. If the set of all {X^ } r >i, is bounded then we do nothing. 
Otherwise we extract from the sequence (X r ) r >i a subsequence (X' r ) r >i, which 
is unbounded on the coordinate i\ also and such that, for every r > 1, X£ < 
X'^ +1 . By applying this argument finitely many times, we can obtain a subset 
/' of I and a sequence a = (Y r ) r >i of solutions such that: 

• for every i £ I', the set {cf} r >i is unbounded; 

• for every r and for every i £ I', Y[ < Y[ +1 \ 

• for every i £ J' — {1, . . . , n} — I', the set {of } r >i is bounded; 

Therefore we can extract from the sequence a a subsequence {Y r ) r >\ such that, 
for every coordinate i £ J', the sequence (Yf) r >i is convergent. Denote by \\v\\ 
the Euclidean norm of a vector v. By the latter condition, for every e > 0, we 
can find two positive integers j < j' , such that: 

VzeJ', \\?f-Yi\\<e. (36) 

Let us consider the vector Z such that: 

V i € J', Zi = 0, V i £ J', Z, = 9? - Y(. (37) 

By Eq. (f3"6")l , one can easily check that ||AZ|| < ce for some fixed constant 
c. Observe now that since a is unbounded on each coordinate i £ V we can 
assume that Zi > 1, for every i £ I' . Let D = J2i'ei' ^'i- ^ ^ s eas ily seen 
that, by replacing in Eq. ([37)) . each Z^ i £ J', with Zi/D, we get a new vector 
Z such that the inequality ||AZ|| < Dee still holds. Since AZ £ Ha and e is 
arbitrary, this contradicts Corollary 0J This ends the proof. 

Now we want to determine a suitable hypercube which includes the set <SU,b- 
To this aim we consider the following linear programming that we call problem 
(Max A, b): 

• maximize Xj, where i is a given index, subject to: 

• AX — b, and 

• X > 0". 
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We make use of the following theorem (cf [TJ], Theorem 2.2) 

Theorem 7 Let B be a matrix in *Z txn , with n > t and let f,d be integer 
vectors. Consider the linear programming problem in standard form given by: 

• minimize d ■ X = X)™=i diXi, subject to: 

• BX = /, and 

• X > 0". 

If an optimal solution exists then there exists also an optimal solution X such 
that for every i, with i — 1, ...,n, \Xj\ < h(3, where h is a non-negative in- 
teger constant, depending only on the matrix B and the vector d, and (3 = 
max{|/j|, \di\}. 

We can now apply the previous theorem to solve Problem (Max A, b). 

Corollary 5 Let A be a matrix in Z txn , with n>t, that satisfies < f53|) . If b is 

a vector in 1} , then there exists a non-negative integer constant Ka, depending 
only on the matrix A, such that the set SA,b *s contained in the hypercube Xi < 
IiaP, for i = 1, ...,n, where (3 = max{|fej|}, with j = 1, 

Proof. By Proposition [U the set Sa, b is bounded in so that Problem 
(Max A, b) admits an optimal solution. Then the claim follows by applying 
Theorem [7] with B — A, f = b and choosing any vector d such that d € {0, —1}™ 
with exactly one non-null component . □ 

Now we give a new proof of Theorem [5j For this purpose, we need some 
preliminary lemmas. The system AX = b is written in the form: 

ana;i + a\2X2 H + a ln x n = bi 

a 2 ixi + a 2 2%2 H + a,2 n x n = b 2 

; (38) 

a t \Xi + a t 2X2 H + a tn x n = b t . 

Let A : 1} — > N be the map defined as: 

Vb = (&i,..., b t ) GZ*, A{b) =h A max 

\<i<t 

where Ha is the positive integer defined in Corollary [5l Consider the family II 
of (hyper-)planes defined by the following equations: 

• for every i = 1, . . . , t, xi = 0, 

• for every i, j with I <i < j <t, Xi — Xj = and Xi + Xj = 0, 
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and let 7Z be the family of regions of 1} defined by II. The first two lemmas 
easily derive from the definition of TZ. 

Lemma 15 For any region R of TZ, there exist a\, . . . ,att G an d j with 

1 < j < t, such that, for every b = (bi, . . . , bt) G R: 

• Vi = l,...,t, \bi\ = (-l) ai h, 

. A(6) = WfhAbj. 

Remark 6 Lemma IT51 shows the role of the family II. Indeed, if b = (&i, . . . , bt) 
is a vector of Z , for every components bi,bj of b, one can determine whether 
\bi\ = bi or \bi\ — —hi and whether < \bj \ or \bi\ > \bj\. 

Lemma 16 Let R be a region ofTZ. There exist a±, . . . , a t S {il} such that, 
for any b = (b\, . . . ,bt) G R, the system i38\) is equivalent to the following: 

\ a n x i + a i2 x 2 H \-ai n x n \= \bi\ 

|o 2 iXi + a 22 X2 H ha 2 «^«|= I&2I 



|a t ixi + a t2 X2 H h a tf! ,a: n | = |6 t |. 

(anxi + 012X2 H h ai„x„) > 

(-1)" 2 (a 2 ia;i + a 22 a;2 H h a 2n x„) > 

< ; ; (39) 

(-l) at (a t ixi + a t2 x 2 H h a 4 „x„) > 

xi > 
x 2 > 



x„ > 
v — 

Lemma 17 Lei R be a region of 1Z. There exists a semilinear set Z of Z* +rl 
such that, for every b = (b%,...,bt) G R: 

C A (b) = G(\b 1 \,...,\b t \,A(b),...,A(b)), 

where G = Gz, t lt t 2 * s the generalized growth function of Z with t\ — t and 
t 2 = n. 

Proof. Let R be a region of 1Z. By Lemma [T51 there exist a%, G {±1} 

such that, for any b = (pi, . . . , b t ) G R, the system (|38[) is equivalent to a system 
of the form (IST))) . Let Z be the subset of all elements of Z t+n defined as: 

Z = {(zi, . . . , z t +n) ■ Zi > 0, i = 1, . . . , t + n}, where 
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• V i = 1, . . . Zi = (-l) ai (ajixi + a l2 x 2 H h a in x„), 

• Vi = t + 1, ...,t + n, Zj = Xi_ t . 

It is easily seen that Z is semilinear in Z t+n . Let G = Qz, t ± ,t 2 be the gen- 
eralized growth function of Z with t\ = t and £2 = n. From the defini- 
tion of Z and by applying Corollary [5] to (|38|) . for every 6 £ i?, one has 
C A (6) = G(|6i|,...,N,A(6),...,A(6)). □ 

Lemma 18 Let R be a region oflZ. The function C A is a box spline on R. 

Proof. Let R be a region of 1Z. By Lemma 1161 there exists an index j , with 
1 < jo < t such that for every b — (bi, . . . , b t ) £ R, A(b) — h A \bj \. 

By Lemma [171 there exists a semilinear set Z of Z t+ra such that, for every 
b=(b u ...,b t ) GR: 

C A (b) = G(\h\, \b t \,A(b), . . . , A(6)) = G(|6i|, . . . , |6t|,^|^„|, ■ ■ ■ , h A \b h \), 

(40) 

where G is the generalized growth function of Z with t\ = t and t% = n. By 
applying Theorem [4] to G, one has that G is a box spline in N* + ™. 

By Lemma O there exists a box spline G' in N* such that, for every b = 

(b u ...,b t )eR-. 

G'(\h\, ...,\b t \) = G(\h\, \b t \,h A \b j0 \, h A \b j0 \). (41) 

We want to show that there exists a F : Z* — ► Z, which is a box spline in 
Z*, such that for every b = (pi, . . . , b t ) e R the following equality holds: 

i^(6 1) ...,6 t ) = G / (|6 1 |,...,|6 t |). (42) 

First of all, recall that by lemma [T5l there exist oti, . . . ,ctt £ {±1}, such that 
for every b = (61, . . . , bt) £ R: 

V^l,...,t, \bi\ = (-l) a *bi. (43) 

Let LTg" = {tti, . . . , 7r m } be the family of planes of K* associated with G' . 
Recall that He satisfies the following properties: 

• TIqi includes the coordinate planes, that is, the planes defined by the 
equations xe — 0, £ = 0, . . . , t; 

• every plane of lie passes through the origin. 

Now let tt(xi, . . . , x t ) = 53i=i t Pi x i — be a plane in Hqi . Then the plane: 

7r'(xi,...,a; t )=E i = w (-1)^/3^ = 

is a plane of K* through the origin. We define as 11^ the family of all such 
planes. It is obvious that all coordinates planes belong to 11^. 
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Let now x — (x\,...,Xt) be a point in R. Since G is a box-spline, G 
associates to the point x* = ((— l) Ql 5;i, (— l) Q *a; t ) a unique quasi-polynomial 
p(x). 

We want to show that the region of x w.r.t. the planes in 11^ determines 
p(x) univocally. Let ir = J2i=i t Pi x i = be a plane in lie ■ Let n' be the 
corresponding plane of n^. Then it is obvious that tt(x*) > (or n(x*) = 0, or 
n(x*) < 0) if and only if n'(x) > (or, respectively, tt'(x) = 0, or, respectively, 
tt'(x) < 0). Then the position of x w.r.t. 11^ determines p{x) univocally. 

Now, let d > the period of p(x). If a given point x — (x\, ...,Xt) in N', 
gives rise to the remainders (di, . . . , d t ) modulo d, then x* determine the same 
remainders (d\, . . . , d t ) modulo d. 

Therefore we let correspond to p{x) the quasi-polynomial q(x) : Z — > Z, 
which has period d and to each sequence of of remainders (di, . . . , dt) modulo 
d associates the polynomial Q(di,d 2 ,---,<it) ■ ^* — * ^> w ith rational coefficients, 
defined by: 

<l(d u d 2 ,-,d t ) (xi, -,x t ) = P(d 1 ,d x ,-,d t )((-i) ai xi, (-l) at x t ) 

The box-spline F is therefore completely specified and this ends the proof. 

□ 

As a corollary of the previous lemma, we obtain. 

Theorem 8 The function Ca of the Diophantine system H38\) is a box spline in 
Z t+n . 

Proof. By Lemma 1181 the function Ca is a box spline on every region of 1Z. 
From the latter fact, one derives that Ca is a box spline on Z* + ™. □ 

Acknowledgements. Special thanks to Corrado De Concini for very useful 
comments and discussions concerning the results presented in this paper. 
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